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O I Abstract. In order to characterise non-equiUbrium growth processes, we study the 

behaviour of global quantities that depend in a non-trivial way on two different times. 

We discuss the dynamical scaling forms of global correlation and response functions and 
^ ' show that the scaling behaviour of the global response can depend on how the system 

^ . is perturbed. On the one hand we derive exact expressions for systems characterised 

"t^ I by linear Langevin equations (as for example the Edwards- Wilkinson and the noisy 

Mullins- Herring equations) , on the other hand we discuss the influence of non-linearities 

on the scaling behaviour of global quantities by integrating numerically the Kardar- 

'^ I Parisi-Zhang equation. We also discuss global fluctuation-dissipation ratios and how 

pi • to use them for the characterisation of non-equilibrium growth processes. 
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1. Introduction 

Non-equilibrium growth is ubiquitous in nature and is encountered in fields as diverse as 
materials science or biological physics [H 121 13] . Kinetic roughening of interfaces displays 
a high degree of universality which has been the focus of many theoretical studies [H |5]. 
Systems with non-equilibrium growth typically have three different dynamic 
regimes. Starting from a flat substrate, the surface initially grows in an uncorrelated way, 
yielding what is sometimes called the random deposition (RD) regime. As time goes on, 
correlations build up, leading to the correlated regime. Even though the nature of this 
regime depends on the rules of the growth process, one generically finds that the time 
dependent mean interface width W{t) increases as a power of time, W{t) ~ t^, where /3 
is the so-called growth exponent. This correlated growth continues until the saturation 
regime is reached for which ly ~ L" where L is the linear size of the system and a 
is the roughness exponent. The different growth universality classes are characterised 
by different values of the exponents a and /3. In the past most studies focused on the 
scaling behaviour of the surface width, as exemplified by the celebrated Family- Vicsek 
scaling relation [El [7]. It is only rather recently (however, see [8] for an early exception) 
that the study of correlation and response functions has been shown to yield interesting 
insights into non-equilibrium growth processes and the related interface fluctuations 

p[ini [mini El [HITS]. 

Two-time quantities are nowadays routinely studied in the context of non- 
equilibrium relaxation and ageing phenomena [IE] where in most cases the focus is on 
local quantities. For example when studying a magnetic system one usually investigates 
the spin-spin correlation function or the response of a spin to a local magnetic field. 
Incidentally, the study of ageing in magnetic systems has revealed that additional 
insights can be gained by looking at global quantities, as for example the magnetisation- 
magnetisation correlation or the response of the magnetisation to a spatially constant 
magnetic field [HI [HI [HI [201 |2l] • 

In the studies of growth processes mainly local two-time quantities have been 
investigated in the past [H [H [TOl [HI [121 [HI [El [H]- Examples include the two-time 
height-height correlation function [SI [9l [HI [HI [HI [H] or the response of the height to a 
local perturbation [HI [HI [HI [El [H] • Some attention has also been paid to slightly more 
complex quantities as for example the two-time roughness or the two-time incoherent 
scattering function [HI [H] • 

Both the local height-height autocorrelation function Ce{t, s) and the local 
autoresponse function Ri{t,s) display simple ageing scaling forms [H [HI [HI [15] . 
irrespective of whether the system is linear or non-linear: 

Ce{t,s) = s-'^fc,{t/s) , Re{t,s) = s-'-'''fR^{t/s) (1) 

where the scaling functions /q and fji^ are power-laws in the long time limit, i.e. 
fceiy) ~ y~^c,i/z g^j^^ fRiiy) ~ y-^R.t/^ for y ^ 1 [H]- The values of the dynamical 
exponent z as well as of the scaling exponents 6^, a^, Ac,^, and Xr^i depend on the 
dynamical universality class. For growth processes the dynamical exponent z is related 
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to the roughness exponent a and the growth exponent (3 by the relation z = a/ (3. As 
shown in the studies [H El [TOl [15] the local autocorrelation and autoresponse exponents, 
Xc,£ and Xr^£, are identical. In addition, they are given by Xc,e = Xr/ = d in systems 
described by a linear stochastic equation. Here d is the dimensionality of the substrate. 
In addition, the scaling exponent of the correlation function, hi, is given by hi = —2a/z, 
whereas for the scaling exponent of the response function the relation a^ = d/z — 1 can be 
conjectured. As for the linear stochastic equations we have that a = {z — d)/2, it follows 
that in linear systems both scaling exponents have the same value, ae = hi = d/z — 1. 
This is different for the non-linear Kardar-Parisi- Zhang equation [23] where for a one- 
dimensional substrate we have a = 1/2 and z = 3/2, yielding bg = —2/3 and ag = —1/3 

m- 

In this paper we study the ageing behaviour of certain global two-time quantities, 
namely the correlation function of the squared width and the response of the squared 
width to a global perturbation. We show that the scaling behaviour of the global 
response depends in linear systems on how the system is perturbed, yielding different 
results for different protocols. This observation is of interest as the global response 
should be readily accessible in smoothening experiments pH], for example. Exploiting 
the fact that exact solutions can be obtained for growth processes described by linear 
Langevin equations, we comprehensively study the scaling behaviour of these global two- 
time quantities for the Edwards- Wilkinson p^ and the noisy MuUins- Herring equations 
[23], thereby distinguishing between various limiting cases. In order to gain some 
understanding of the behaviour of these quantities in non-linear systems, we also discuss 
some data that have been obtained by numerically integrating the non-linear Kardar- 
Parisi-Zhang equation [24]. Our results indicate that in the non-trivial correlated regime 
the two-time global quantities generically exhibit a behaviour of full ageing. We also 
discuss the global fluctuation-dissipation ratio and show how this quantity can be used 
for the characterisation of growth processes. 

Our paper is organised in the following way. In the next section we remind the 
reader of the typical behaviour of a growing surface. This is done with the help of 
a microscopic deposition model [7] that is very well described by the one-dimensional 
Edwards- Wilkinson equation. This model also allows us to motivate the different global 
perturbations that we are going to discuss in the following sections. Section 3 is devoted 
to the exact computation of the global two-time quantities in cases that are described 
by linear stochastic Langevin equations. The exact results allow us to comprehensively 
investigate all possible dynamical regimes. In Section 4 we present our results for the 
one- dimensional Kardar-Parisi- Zhang equation where we highlight some commonalities 
with and differences to the linear systems discussed in the previous section. Finally, 
Section 5 gives out conclusions. 
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2. Motivation: A deposition model 

In the following we briefly discuss the behaviour of a simple deposition model that turns 
out to be an excellent representative of the Edwards- Wilkinson universality class. This 
model is then used to motivate the different protocols discussed in the next sections for 
measuring the global response. 

2.1. Definition of the model 

Our deposition model [7] is based on Family's original random deposition with surface 
relaxation (RDSR) process [25] and differs from this model by the diffusion step. In 
the RDSR process a particle deposited on the surface is allowed to jump to one of the 
neighbouring sites if this site has a lower height than the site of deposition. In our 
model we assign an energy Ei{t) = g hi{t) to the column at site i where hi{t) is the 
height of that column at time t. The constant g can be thought of as the gravitation 
constant, for example. Starting from an initially flat substrate, particles of mass one are 
deposited on randomly chosen sites and then allowed to diffuse locally after deposition. 
For a diffusion step taking place at time t, we select one of the neighbouring sites j 
at random and accept the jump with the temperature and time dependent (Metropolis 
like) probability 

' 1 if Ej{t) < Eiit) 

\ if E,{t) > E,{t) 

In the following we choose units thus that the Boltzmann constant fc^ = 1. 

In contrast to the original model there is a non-vanishing probability that a 
deposited particle jumps to a neighbouring site with a higher height than the deposition 
site. We assume this jump to be thermally activated and to depend on T (the 
temperature of the substrate). As we discuss in the following, the ratio T/g is a 
parameter that governs the morphology of the growing interface and allows us to study 
the response of the surface to a change in external conditions. 

It is instructive to look at the behaviour of the model in the limits of T — )■ and 
T — 7- oo. At zero temperature no jumps to sites with higher height are allowed, and 
the particle is incorporated into the aggregate at the selected neighbouring site if the 
column at that site is shorter than at the initial site. Thus, for T — t- we recover the 
RDSR process [25]. In the opposite limit, T — ;■ oo, however, a particle will always jump 
to the selected neighbouring column, irrespective of the height difference. As a result, 
the different columns will grow independently, yielding an uncorrelated surface as for 
the RD process. For intermediate temperatures, a crossover between the RD and the 
RDSR processes is observed, with the crossover point depending on the temperature. 

The temperature dependence of the model is illustrated in Figure [H Suppose a 
particle is deposited on top of the middle column in the conflguration shown, the plot 
provides the T/g dependence of the probabilities for having this particle end up at one 
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of the three sites. In the original Family model the particle would always come to rest 
on top of the left column. 




T/g 

Figure 1. Probabilities that for the shown configuration the particle, initially 
deposited on top of the middle column, comes to rest on top of one of the three 
columns. 

Of course, temperature has been introduced in theoretical studies of growth 
processes prior to our work (see [26], [271 EHl EHl [30] for some examples). For a typical 
microscopic model studied in this context one often assumes that the atom's hopping 
follows an Arrhenius-like rate that is proportional to e"^"/'^^"^, where the activation 
energy En is itself proportional to the number of bonds formed by the atom before 
the hopping attempt. Whereas this is surely a rather realistic modelling of diffusion 
processes on a crystal surface, we have opted here for a simpler approach where the 
probability for a particle to hop depends on the height difference between the actual 
site and the proposed new site. (Note that jumps to sites with higher heights have 
also been allowed in other microscopic growth models [261 IS])- As we will show in 
the following, all aspects of our simple model agree perfectly with the solution of the 
stochastic Edwards- Wilkinson (EW) equation with a single fit parameter. 



2.2. Interface width 

In the following we are interested in the dependence of the surface width W on time t, on 
the ratio T/g, and on system size L (in this section we only discuss the one-dimensional 
case). In our simulations, we have a wide range of T/g values and several L's up to 
1000. For simplicity, we always start from a fiat surface, i.e. W = ioi t = 0. For the 
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data discussed below, the unit of t is one Monte Carlo Step (i.e., L particles deposited) 
and we averaged over 1000 independent runs with different random numbers. 

In Figure [2t,c we show the time dependence of the surface width for, respectively, 
the case with fixed L = 1000 at different T/g^s and the case with a fixed T/g = 1 and 
various L's. As for the RDSR process one distinguishes three regimes separated by two 
crossover points: a random deposition (RD) regime, followed by a EW regime, with 
a final crossover to the saturation regime. In contrast to Family's original model, the 
initial RD process is not confined to very early times t < 1 but extends to larger times. 
In fact, the crossover time ti between the RD and the EW regimes is shifted to higher 
values for increasing temperatures and diverges in the limit of infinite temperatures. 
As the crossover is smeared out, we identify the crossover point with the intersection 
point of the straight lines fitted to the two linear regimes in the log-log plots. Due to 
the nature of the uncorrelated deposition of particles (Poisson process), we have the 
identity W"^ = t in the RD regime, yielding the width Wi = ^Jt[ at the crossover point. 
In the EW regime the relation between width and deposition time changes to VT oc t^^^. 
This regime extends up to a second crossover point (^2, ^^2)? whose precise location also 
depends on the value of T and beyond which the final saturation regime prevails. The 
crossover between the different regimes is further illustrated in Fig. WpA where we show 
the time evolution of the effective exponent 

for the two cases. 

Next, we turn our attention to the stochastic EW equation (see also the next 
section) 

^^ = z.V^Mx,t)+r^(x,t) (4) 

where //(x, t) is a Gaussian white noise with zero mean and covariance (^^(x, t)ri{'y, s)) = 
Z)(5'^(x — y)S{t — s) {d is the dimensionality of the substrate) and z/ is the surface tension 
or diffusion constant. This equation can be solved exactly to give us in one dimension 
the width (squared) 

where g„ = 27m/ L and the sum is over [— L/2, L/2] but excluding the zero mode: n = 0. 
The effective exponent is therefore given by 



utJ2e' 



-2vtq, 



^•^ff E [1 -e-^^^'^l] /ql • ^^^ 

n 

In the expression ([5]), we must fix Z^ = 1 in order to agree with W'^ = t in the RD 
regime. This leaves us with just one free parameter, namely, v{T/g)^ which can be 
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Figure 2. (a) Log-log plot of the surface width vs time for a system of size L — 1000 
and different values oiT/g. The dashed lines have the slopes 1/2 and 1/4 expected in 
the random deposition and EW regimes, respectively. The location of both crossover 
points depend on temperature. The full lines are obtained from fitting the exact 
solution of the EW stochastic equation. Here and in the following error bars are 
smaller than the symbol sizes, (b) Time evolution of the effective exponent (|3]) for the 
data shown in (a), (c) Log- log plot of the surface width vs time for systems of different 
sizes evolving at the value T/g = 1. (d) Time evolution of the effective exponent (j3|) 
for the data shown in (c) . The full lines are derived from the fits to the exact solution 
of the EW stochastic equation. 



obtained by fitting the numerical data to the theoretical expression, see Table [H The 
result of this procedure can be summarized by 

with a = 4.23 and b = 2.13. As the solid lines in Fig. [2] show, excellent fits are achieved 
with these values of u. Clearly, the theoretical curves are in very good agreement with 
the data over the entire range of L's and T/g^s explored. Note in addition that it is 
possible to collapse all data onto a single curve [7]. 

It is worth mentioning that a 1/T dependence of the surface tension u is 
typically encountered in experiments on step fluctuations [32l |33l EH [35] or on surface 
smoothening [Tl] which are described theoretically by linear Langevin equations. 

The evolution of the width shown in Fig. [2] is very reminiscent of the behaviour 
encountered in competitive growth models [36l|371[3Hll39lll0lllIlll2lll3lllIlll5lll6]. 
In both cases two different scaling regimes are found whose ranges depend on a system 



Characterisation of non- equilibrium growth through global two-time quantities 



T/g 


0.1 


1 


10 


100 


V 


0.2300 


0.1768 


0.0379 


0.0046 



Table 1. Values for the surface tension v at different values of T j g that result from a 
fit of the numerical data to the exact solution of the Edwards- Wilkinson equation. 



parameter. In the competitive growth models one considers a mixture of two different 
deposition processes where one of them takes place with probability p whereas the other 
takes place with probability 1 — p. One example is the RD/RDSR model [36] where 
the deposition happens according to the RDSR rules with probability p and to the RD 
rules with probability 1 — p. Whereas ioi p = 1 and p = only one of the processes 
is realized, for general values of p the mixture of the two processes leads to a crossover 
between the two regimes where the crossover time and width depend on the value of p. 
In our model, the ratio T/g plays the role of the quantity p, the main difference being 
that we do not artificially choose between the different processes, as in our system the 
competition is intrinsic and governed by the value of the temperature. 



2.3. Global perturbations of the growth process 

Experimentally, a quantity that can be changed easily and that gives way to a global 
perturbation is the temperature. If an experimental system is described by a stochastic 
Langevin equation, temperature can enter either through the surface tension u (also 
called diffusion constant or mobility, depending on the physical context) or through 
the noise. In the next sections we shall study the responses to two different global 
perturbations: either we keep the noise unchanged and suddenly change z/ (as it is 
the case in surface smoothening experiments or in step fluctuation studies) or we keep 
the surface tension constant and change the noise (as it is the case in some deposition 
experiments where the noise in the particle flux can be changed experimentally). Both 
protocols have been used for studying the change in morphology of a growing surface 
when changing experimental conditions ^7\ [13]. 

We can use our deposition model in order to see which global quantity is best suited 
to study the changes due to such a global perturbation. Let us consider at temperature 
T a configuration where the height of the deposed column at site i is hi. Setting the 
lattice constant in vertical direction to be 1 and setting the potential energy of the initial 
fiat surface to be zero, the potential energy f/j stored in a column of height hi is (with 
the mass of a deposed particle set to 1) 



^ = 1 + 2- + 

rji rji ' rji ' 



+ {h-l)^ = hik-l)^. 



(8) 



N 



Shifting the value zero of the potential energy to the average height h = jf J2 hi, where 



i=l 



A^ is the number of sites on the substrate, we obtain for the total potential energy of 
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our configuration the value 

^ = E I = 1= E {h<ii'. - 1) - MS - D) = |; E {1'^ - h') = f i^^P) 

2=1 i = l 2 = 1 

where W'^ = jqjl [hf — h ) = j^ J2 [hi — hj is the squared width. From this it follows 

i=l i=l 

that yiy^ is the quantity conjugated to g/T and, as g/T ~ i/, to z/. This motivates the 
use of the square of the surface width (instead of the surface width itself) as the global 
quantity to be studied in the following. 

3. Linear Langevin equations 

Linear Langevin equations jl] are discussed in a variety of physical situations, ranging 
from equilibrium step fluctuations [Ml SS |35l [31] to film growth [501 EB [52] and from 
magnetic systems ^16] to elastic lines in a random environment |53]. In the following we 
focus on two simple but generic cases which can be summarised by the equation 

^h^ = -u{zVrh{^,t) + r/(x,t), (10) 

where m is an even number and the noise is uncorrelated and with zero mean, 

(r/(x,t)r^(x',t')) = ^'5'(x-x')5(t-t') . (11) 

For m = 2 we recover the well-known Edwards- Wilkinson (EW) equation ^22j and for 
m = 4 we obtain the noisy Mullins- Herring (MH) equation [23]. Note that for these 
systems the dynamical exponent z = m. 

As already discussed in the introduction, some aspects of ageing in these systems 
have been elucidated previously through the study of local quantities. In the following 
we are going to characterise this relaxation through global two-times quantities related 
to the square of the surface width. 

3.1. Global two-time correlation function 

The first global quantity that we are going to discuss is the connected correlation 
function (here and in the following (■ ■ ■) indicates an average over the noise) 

C{t, s) = {W\t) W\s)) - {W\t)){W\s)) (12) 

where 

W'{t) = ^W\t) = \'Z(h{^,t)-hf (13) 

W\t) = ^^Y.{K^^t)-hf (14) 

is the square of the surface width at time t on top of a d-dimensional substrate of linear 
extension L. Here x labels the different substrate sites. We assume in the following that 
t > s and call s the waiting time and t the observation time. 



and 
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Writing both the height /i(x, t) and the noise 77 (x, t) as a sum over reciprocal lattice 
vectors, 

^(x,t) = X!^qWexp(iq-x) , r7(x,t) =^r7q(t)exp(iq-x) , (15) 

q q 

equations (TTOj) and ( ITTj) become 

^ = -i^cTh^ + v^ (16) 

and 

{UtW^it')) = DL-X^^'S{t-t'). (17) 

The general solution of ( TT6l) is readily obtained to be 

/iq(t) = exp(-z/q'"t)/iq(0) + rrfrexp(-z/q"(t-r))r/q(r). (18) 

^ u 

Assuming an initially flat surface, /iq(0) = for every q, one obtains the following 

two-point correlation function of the surface height in q-space: 

{Kit.)Kih)) = ^--'^'''''^'^'''-^ (e^(c..r)*< _ 1) si,^^, (19) 

where g^ = |qj| and t< is the smaller of the times ti and ^2- This well-known result will 
be useful in the following. 

For the computation of the correlation function (TT2|) we remark that the square of 
the surface width ( !T4l) can be written in terms of the Fourier components of the height 

as 

^'W = ^E(EM^)exp(zq-x)] (20) 

^ X Wo / 

where we took into account that the average height of the surface is equal to /iq=o- This 

then yields 

{W\t)W'{s)) = i 5: Y: {K{t)h^'{t)h^{s)h^is)) 

x,x' q,q',p,p'5!^0 

exp [i(q + q') • x] exp [i{p + p') • x']. (21) 

The calculation of (TT2|) is therefore reduced to the determination of the four-point 
correlation function of surface height in q-space. Using the Gaussian statistics of the 
linear equations, this four-point function can be expressed through two-point functions 
of the form ( fT9l) , yielding 

(V(^l)^q2(^2)V(^3)V(^4)) = (V(^l)^q2(^2))(/iq3(i3)V(^4)) 

+ (V(^l)V(^3))(/^q.(t2)V(i4)) 

+ (V(il)V(^4))(/iq.(t2)/iq3(t3)). (22) 

Inserting this into ( 1211) and using the identity 



'^ '27dnm\ 



n=0 



E e^P [^^) = NSrrr,o (23) 
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yields the expression 

{W\t)W\s)) = ^ I E -^— (l - e-'^'"'^*) (l - e-'^"^^' 

1 

-,2r?i ' 



I y^ ^ ~2q'^ut t 2q^us , ^2uq^ 



The first term on the right-hand side being just (iy^(t))(W^^(s)), we finally obtain 

C«. -) = S E i^-^'"-' (^''"" + <''^"" - 2) ■ (24) 

The behaviour of C(t, s) is controlled by two length scales: It = i^vt^^^ and 
/s = (2z/s)^'"^, with /( > Ig. Depending on the values of It and Is and their relations to 
the maximum and minimum values of g, qmax = vrv5 and qmm = 2tt/L, different regimes 
are encountered. For example, if at time t the length It < l/qmaM^ one still is in the 
short time regime where the random deposition process prevails. On the other hand, 
if at time t one has that k > l/qmin, one is the saturation regime. Finally, the system 
is in the correlated regime when l/qmax < U < ^llmin- We discuss in the following the 
functional form of the correlation function for the different regimes (see Table [2] for a 
summary) . 

1. If It < ^jqraax wc are in the short time regime. In that case all exponentials in flM|) 
can be expanded and only the leading terms need to be retained, i.e. e"^''™'^* ?« 1 — 2g™z/t 
and e^i"^^^ + g-2ij"Vs _ 2 ^ (2g'"z/s) for all q. Replacing the sum by an integral, we 
obtain 

n2o2 / T \d T, 

C{t, s)^^ (- j n, J^ q'-' (1 - 2q-ut) dq 

D^TT'^/^L'' of 2d „, \ 
^ ,^, ,, , , s^ 1 - 7r'"z/t , 25 

where ild = 2'k'^/'^ /V{d/2) is the solid angle of the (i- dimensional sphere. Therefore, in 
the short time regime the correlation function decreases linearly with time. 

2. If the final time is in the correlated regime, i.e. if l/qmax < k < l/qmin, one 
has to distinguish between two different cases, depending on whether Ig < l/qmax (the 
system was in the RD regime at time s) or l/qmax < h < k (the system was in the 
correlated regime at time s). In the first case we can again replace e^'^™'^'^ + e"^''™'"* — 2 
by {2q"^vs) and obtain the expression 

= A„,4'H-''I'\ (26) 

X By It < l/qmax we mean that the order of It has to be much smaher than the order of l/qmax- There 
are comphcated crossover effects showing up when It and l/qmax are of comparable magnitude that we 
are not going to discuss. Of course, the same remark appHes when comparing one of the length scales 

to l/qmin- 
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where 

^ D^T{d/m)L^ 

yielding a power-law decay with an exponent d/m. In the second case, which is the most 
interesting one as both times are in the correlated regime, we expand all exponentials 
in fl2^ and integrate term by term. This yields 

C(t, s) = A^,,s' \ t~^i'^ + ^^^t-"'"^ + T^^t--^/"^ + ■ ■ -I 




^ Am^sh-"'"' (28) 

where the last line gives the asymptotic form for t ^ s. Interestingly, we recover 
asymptotically the same leading behaviour as for the case Is < l/qmax, i-e. in the 
long time limit any memory of the regime that prevailed at time s is lost. Looking 
closer at the expressions ( l26l) and ( l28l) . we see that in both cases the global correlation 
function can be cast in the standard form |l6j C{t, s) = s~^fc{t/s) where b = d/m — 2 
and fc(t/s) is a scaling function that only depends on the ratio t/s. In addition, the 
scaling function is asymptotically given by a power-law, fcit/s) ~ {t/ s)~^^^^ ^ with the 
autocorrelation exponent \c = d and the dynamical exponent z = m. In Figure [3H,b we 
compare the correlation function (IMl) . where the sum has been evaluated numerically, 
with the asymptotic power- law fl28|) . see Figure |3^ for the EW case and Figure [3]d for 
the MH case. In all cases the asymptotic regime is accessed very rapidly. 

We can now also compare for this regime the global correlation function with the 
local height-height correlation, see [HI [9l [10]. For the local correlation function one 
obtains the scaling exponents bi = d/m — 1 and Xc/ = d. Thus, for both the local and 
global quantities the long-time behaviour is governed by the same power-law exponent, 
Xc = Xc,e, but the value of the scaling exponent b is different from the value of bi. This 
yields interesting differences at the upper critical dimension which is d = 2 for EW 
and d = 4 for MH. Indeed, whereas the local exponent bi then vanishes and the scaling 
function of the local correlation shows a logarithmic dependence on time (which makes 
the extraction of the correct scaling behaviour from experimental or simulation data 
rather tedious) [9], the global exponent remains different from zero and no logarithms 
are occulting the scaling behaviour of the global correlation, making this a much better 
quantity for studying dynamical scaling close to the upper critical dimension. 

3. In case the observation time is in the saturation regime with It > 1/qmin, the 
leading contribution is given by the term in e"^''™™'^*. If the waiting time is also in the 
saturation regime, we have Is > l/qmin and therefore 

dD^L'^"' 
4(27r)2™z/2' 



Cit^s)^ ,,, ,,^ e-^C..K^-^), (29) 
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SR 


— 


— 
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Table 2. Summary of the asymptotical dependence of the global correlation function 
for the different cases. RD corresponds to the random deposition regime, CR to 
the correlated regime, and SR to the saturation regime. The values of the different 
coefficients can be read off from the equations in the main text. 



i.e. the global correlation is time-translation invariant. This is of course expected, as 
the saturation regime corresponds to the steady state regime. If the waiting time is 
such that Is < l/qmin, then we can approximate for q = q^in the term involving s by 
(2g™^„z/s) which yields 

C{t, s) ^ c/D2^2g-2C,„^'t _ ^3Q^ 

It is worth noting that in both cases the decay time in the exponential, r = l/(2g™j„z/), 
depends on the diffusion mechanism through the values of u and m but does not depend 
on the dimensionality of the substrate. 



3.2. Global two-time response functions 

In order to measure the global response we are going to perturb the growing system 
either through a change in the surface tension i/ or in the strength D of the noise 
correlator, as discussed in Section 2.3. Having thus perturbed our system, we then 
study its relaxation to the state that is realised when the system evolves all the time at 
fixed values of u and D. 

In the first protocol, we start from an initially flat surface and let the system evolve 
with a surface tension /i until the waiting time s after which we set the surface tension 
to the final value u. In order to monitor the relaxation of the system for t > s, we 
compute the function 



X{t,s) 



{W'),^,{t,s)-{W^).{t) 



(31) 



where the index " fi ^ u" indicates the perturbed system, whereas the index "z/" stands 
for the system that is evolving at a fixed value u of the surface tension. The quantity 
e = u — n measures the strength of the perturbation. Usually one is interested in small 
perturbations for which |e| ^ 1, and we will start with this in the following. However, 
as we discuss later, the response (13T]) displays simple ageing and dynamical scaling even 
for very strong perturbations. 

This protocol has been used previously in [13] in order to study the asymptotic 
long-time behaviour of a one- dimensional perturbed EW system. No attempt was made 
in that work to determine the scaling functions in the ageing regime. 
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Figure 3. Global correlation (a,b), global response to a change of the surface tension 
(c,d), and global response to a change of the noise strength (e,f) when both times t and 
s are in the correlated regime. The symbols are obtained by numerically evaluating 
the exact expressions (|24l) . psp . and (jH]) whereas the lines indicate the asymptotic 
power-laws (|28|) . (p9|) . and (|48p . Panels (a), (c), and (e) show data for the EW case 
and dimensions d = 1,2 (from top to bottom), panels (b), (d) and (f) show data for 
the MH case and dimensions d = 1, 2, 3, 4 (from top to bottom). The different symbols 
correspond to different waiting times where sq = 1000 below the critical dimension 
and So = 200 at the critical dimension (which is d — 2 for EW and d = 4 for MH). 
The linear extension of the system is L = 2^^ for d = 1, L = 2^ for d = 2, L = 2^ for 
d = 3, and L = 2^ ioi d = 4. 



It has to be noted that the response x(t, s) is a time integrated response that gives 
the reaction of the system to a perturbation lasting until time s. This time integrated 
response is related to the response R due to an instantaneous perturbation by x(t, s) = 

s 

J duR(t,u). Assuming a standard scaling behaviour for R, i.e. R(t,s) = s^^^'^fji(t/s), 



it follows that the time integrated response should scale as x(^5 s) = s~"'f^{t/s) [T6] . 
The expression for x(t, s) is readily obtained by remarking that the solution of the 
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Langevin equation becomes for t > s 

h^{t) = e-'^'''"(*-^)/iq,^(s) + J^ dre-'^^'^^'-^^r^^ir), (32) 

where 

Vm(s)= rrfre-^^'"(^-^)77q(r) (33) 



/o 

is the solution of a surface that, starting from a flat initial state, evolves until time s at 
the constant surface tension /i. Inserting this expression into ( l20ll yields 



(W'),^At,s) = -Y.-;^{-e 



e 



4 ^^q^ i/x // 

+ i(l_e-2'?™K*-))| (34) 

for t > s, which finally gives the expression 

X(*. -) = § E ^«-^'-'"- { ^ (l - ^-^'■■'") - ; (l - e-'"""-) ] (35) 

for our response function. 

For small perturbations with |e| ^ 1 we can proceed as for the correlation function 
in the previous section. As we have the same cases and therefore can use the same 
approximation schemes, we shall only quote the final expressions, see Table [3] for a 
summary. 

1. At early times, where t is still in the RD regime, we have a linear decay, 



^^^'^^'^2'^r(d/2)(d + m)^' 

[ r,d + ^ m 

X < 1 - 2 vr 

d + 2m 



^{t - s) + -{ly + fi)s 



+ 0{vtf (36) 



2. When t is in the correlated regime, we obtain for a perturbation that lasted only 
until some time s in the RD regime the following expression for the response: 

x{t,s)^Bm,ds\t-sr'l'^-\ (37) 

with 

^ DdT{d/m)L'' 

If, however, the perturbation ended inside the correlated regime, we have 

X{t, S) = Bm,dS^{t - S)"™"^ 

^ B,^^ds\t - s)-^-\ (39) 

yielding the same leading behaviour as for a perturbation that already stopped in the 
RD regime. In the correlated regime we can therefore cast the response in the standard 
form xi^is) = s~°- fy^{t/s) with a = — — 1 and fxit/s) ~ (t/s)~'*'«/^ for t/s large, with 
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Table 3. The same as Table [21 but now for the response (|35|) to a smah perturbation 
during which the surface tension is changed. The values of the different coefficients 
can be read off from the equations in the main text. 



Xr = d + m. We show this scahng behaviour in Figure |3] c,d for the EW and for the 
MH case and compare the scahng function to the asymptotic power-law. 

One can now also compare the values of the exponents with those obtained for 
the global correlation function ( TT2l) as well as with those obtained for the response of 
the height to a local perturbation, see [9]. Comparing with the global correlation, one 
notes that we have the interesting situation that a ^ h and Xr 7^ Xc- It is rather 
uncommon that the response of some quantity and the correlation of the same quantity 
display a different asymptotic behaviour in the ageing regime [16] . Mathematically, 
this intriguing observation can be traced back to the fact that in the expression ( 12^ 
for the correlation g^™ shows up in the denominator whereas q^ is encountered in the 
denominator of the corresponding expression ( !35l) for the response function. In order 
to compare the global response with the local response, we need to look at the time 
integrated local response. Inserting the expressions given in [9] for the local response R^ 

s 

into the integral XfX^-, -s) = / duRi{t,u) yields for the local time integrated response the 



form Xiitj ■s) = s'""' fxti'^^ ^) where the exponent a^ has the value ai = ^ — 1, whereas 
the exponent governing the scaling function in the long time regime is Xr^i = d. We 
therefore have that the scaling exponent a = a^ is the same for both responses but that 
the global and the local autoresponses differ, Xr = Xr^ + m. 

3. Finally, if the observation time is in the saturation regime, we obtain that 

Xit,s)^ __^^^_ e-^C.„K^-^) (40) 



2m-d+2^m^j^ 



for a perturbation that lasts until the saturation regime and 



Xit,s) 



I'm 



-s'e 



-2q"^. i^(t-s) 



(41) 



for a perturbation that comes to an end before entering the saturation regime. 

Let us close the discussion of the first protocol by briefly looking at the case of 
strong perturbations. As /i and v then strongly differ, we need to take into account the 
existence of a third length scale /a = (2/is)^/™. This makes the analysis more involved. 
We first note that for the case z/ :» /i we always have the relation It > Ig > l^. It is 
then readily verified that we have basically the same three cases as before, with the 
response given by the expression ( l36l) respectively ( l39ll when the observation time t is in 
the RD regime respectively in the correlated regime. In addition, a simple exponential 
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decay is encountered when t is in the saturation regime. If /i ^ z/, the situation is 
more comphcated, as we do no longer have the above hierarchy of the length scales, 
and in many cases no simple dynamical scaling is observed (see also Fig. 5 in [T3]). 
Of particular interest is the case I3 > l/qmin and l/qmax < h < ^llmin for which the 
observation time is in the correlated regime, as this yields for the global response the 
expression 



X(t, s) = ^Y. e-^^"^^(*-^)- [-us + 2q-{usf + ■ 



This is a remarkable result, as we recover for an arbitrary large perturbation a power- 
law relaxation and simple dynamical scaling. In addition, the exponents describing the 
scaling behaviour have values that differ from those obtained for a small perturbation. 
Indeed, from the expression fH2l) we obtain a = ^ — 1 and \r = d. 

For the second protocol, we proceed analogously. We again start from an initially 
flat surface and let the system evolve with the strength D' of the noise correlation until 
the waiting time s at which we set this strength to the final value D. We then compute 
the function 

e 
with e = D' — D. This is again a time integrated response. A straightforward calculation 
yields the expression 

X«(*.-) = ^i:^«-^'"'" («='"""-!) (44) 

for this global response. Again, we briefly discuss the limiting cases in the following (see 
Table H. 

1. At early times with It < l/qmax we have again a linear decrease of the global 
response: 

2. In the correlated regime we obtain for a perturbation ending before the end of 
the RD regime: 

XD{t,s)^Cm,dSt-''/'^ (46) 

with 

fy _ ^m^ (A^j) 

If the perturbation only ends in the correlated regime, the result is 

XDit, S) = Cm^St-'/"- (^ + i? + ■ ■ ■) ' ^^^) 
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Table 4. The same as Table |3l but now for the response (j43|) to a change in the 
strength of the noise correlation. The values of the different coefficients can be read 
off from the equations in the main text. 



yielding the same asymptotic behaviour. This scahng behaviour is shown in Fig. [3^,f 
for the EW and MH cases. Comparing the expressions (HHll and ( l39l) we see that the 
scahng properties of the response of the square of the surface width depends on how 
the growing interface has been perturbed. More precisely we find that (1) the scaling 
exponent a is independent of the protocol, but that (2) the autoresponse exponent \r 
is now equal to d, i.e. the response of the system decays in the long time limit slower 
when the perturbation is due to a change of the noise strength. Interestingly, the global 
response to a change in D has the same power-law as the global autocorrelation. 

3. Finally, if the perturbation continues until the saturation regime is reached, then 



XD{t,s) 



dU 



.-2CinKt-^) 



(49) 



2™+%™!^ 

If, on the other hand, we measure in the saturation regime the response to a perturbation 
that ended before reaching that final regime, we have that 



XD{t,s) ^ dse ^'i'™-"''* . 



(50) 



3.3. Global fluctuation- dissipation ratios 



The fiuctuation-dissipation ratio pi] has been introduced as a generalisation of the 
celebrated fiuctuation-dissipation theorem to non-equilibrium systems. Whereas this 
ratio has some appealing features, as for example its simplicity and the possibility to 
assign an effective temperature to a non-equilibrium system [55] , many studies in model 
systems have revealed major problems with that quantity, ranging from observable 
dependent effective temperatures [561 EZ] to the appearance of negative temperatures 
[581 [59]. As a result, the usefulness of the fluctuation-dissipation ratio is in general very 
restricted, even so it might occasionally yield interesting insights in speciflc systems. 

The fluctuation-dissipation ratios of local quantities and the related effective 
temperatures have also been studied for the one-dimensional EW equation [IT] (higher 
dimensional cases and the MH cases can be inferred from the equations given in [9])- This 
study shows that the ratio of the response Ri(t, s) of the height to the height-height 
correlation Ce{t, s) allows in the linear systems to characterise the different growth 
regimes through an effective temperature. 
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For a time-integrated linear response x{ty s) ^i^nd the conjugate correlation function 
it is convenient to discuss the ratio 

.m.)^^ (51) 

and measure the limit value X^o = lim^ s-oohnit ^oo X{t,s) [IE]. Strictly speaking, 

taking this limit does not yield non-trivial results for our global quantities. Assuming 
this system to be infinite (such that one remains in the correlated regime in the long- 
time limit), one immediately sees from the Tables [2] and |3] that the trivial value X^ = 
is obtained for the global fluctuation-dissipation ratio, and it is not possible to introduce 
a meaningful effective temperature. Still, the quantity X(t,s) does yield at finite times 
a behaviour that is characteristic for the different regimes, as shown in Fig. H] for the 
ratio formed by the global response to a change in u and the global correlation. In order 
to understand this figure, let us look at a case (the case s = 1 in the figure) where 
the waiting time is in the RD regime, li t ^ s, i.e. t is also in that regime, the ratio 
(15T|) is constant: X = jjf^r^- If t is such that 1/qmax < k < ^/qmin, i-e. t is in the 
correlated regime, X is inversely proportional to t: X = 2Dvm ^~^' '^^^ same behaviour, 
albeit with a different pre-factor, is obtained if we consider the global response to a 
change in D instead. Finally, in the saturation regime, with It > 1/qmin, X is again 
constant, with X = ^ dDL"^"" ■ These three regimes are separated by crossover regimes 
when It ~ l/^max or U ~ 1/q'mm- If we choose a different waiting time, s = 100 or 
s = 10000 in Fig. HJ the ratio X{t, s), after some non-universal behaviour for t — S'^s, 
rapidly evolves towards the same master curve as the s = 1 data. 

4. Non-linear growth 

In order to assess the importance of non-linearities, we also studied global quantities 
derived from the one-dimensional KPZ equation. We thereby complement the studies 
[121 US] which focused on local quantities. 

The KPZ equation has been shown to describe faithfully kinetic roughening and to 
be a paradigmatic model for the description of a large range of non-equilibrium systems 
[2]. This well-known equation is given by the expression 

^^^ = uV'h{^, t) + A(V/i(x, t)f + r^(x, t) (52) 

at 

and differs from the EW equation by the non-linear term proportional to the parameter 

A. Here //(x, t) is the usual Gaussian white noise with zero mean and covariance 

(?7(x, t)r7(y, s)) = D6'^{:s. — y)S{t — s). The numerical integration of (!52|) has been 

the subject of many studies and different discretisation schemes have been proposed in 

order to handle the non-linearity (see, for example, [601 Ell Ell ESI El])- 

We also studied for this system the global correlation and response functions 

introduced in the previous sections. In order to make sure that our results are 

independent of the integration schemes, we computed these quantities using different 

approaches, namely the scheme of Lam and Shin [61] as well as the strong coupling 
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Figure 4. The global fluctuation-dissipation ratio for the one-dimensional EW case, 
with m = I. The ratio ((5T|) displays a different behaviour in the different regimes. 
The data have been obtained by numerically evaluating the exact expressions derived 
in the previous subsections. The parameters are L = 500, D = 1, and ly — 0.001. 



scheme proposed by Newman and co-workers [60l [65] . We carefully checked that we are 
in the correlated regime at both times s and t. We also verified that both approaches 
yield the same exponents and scaling functions for the studied quantities |§|. The data 
obtained from both schemes differ by a non-universal pre-factor, as expected. Small 
deviations between the data sets are observed for t ~ s, but this is again expected as 
for the Lam/Shin scheme we use finite values of A, whereas the Newman scheme works 
in the strong-coupling limit where A — > cxd. 

Fig. [5] summarizes our main findings regarding the scaling behaviour of global 
quantities in the one-dimensional KPZ universality class. For the autocorrelation (iT2ll 
we find the scaling form 



C{t,s) 



fcit/s) 



(53) 



with fdv) ~ y ^^^ for y ^ I. For the response functions we note, and this is in strong 
contrast to the results we obtained for the linear equations, that the scaling behaviour 

§ Due to the very nature of the Newman algorithm, the surface tension is not a parameter that can be 
changed in that scheme. Therefore, we only computed with that method the global correlation function 
as well as the global response to a change in the noise strength. 
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Figure 5. Log- log plots of global quantities as a function of t/s obtained from the one- 
dimensional KPZ equation in the correlated regime: (a) global correlation (b), global 
response to a change of the surface tension ly, and (c) global response to a change of 
the noise strength D. The data in (a) and (c) have been obtained with the Newman 
algorithm [601 ESj , whereas the data in (b) result from the algorithm of Lam and Shin. 
We checked that both algorithms produce the same exponents and scaling functions 
for a given quantity. The system size is L = 10000, with X = 1, ly = I, and D = I. 
For the response shown in (b), v has the value of 1.1 until time s at which point it is 
changed to the value 1, whereas for (c) D was changed from the value 1.1 to 1 at time 
s. The data shown have been obtained after averaging over 5 x 10^ realisations for the 
correlation and 5 x 10^ realisations for the response. The lines indicate the asymptotic 
power-law behaviour. 



does not depend on the protocol that we use for perturbing the system [||J. Whether we 
change the surface tension or the noise strength, we always end up with the following 
scaling behaviour of the global response: 



X{t,s) 



,2/3 



fxit/^) 



(54) 



with fxiv) ~ y^^^^ for y ^ 1. This difference between the global response in the 
linear and the non-linear cases is remarkable. At this stage, we can only speculate 
on the origin of this, but it seems that the restoring mechanism, that is responsible 
for the relaxation of the perturbed but still growing interface, originates in the non- 
linearity, thereby yielding the same scaling behaviour independent of the nature of the 
perturbation. Further studies seem to be in order to fully clarify this point. 

II Noting that by rescaling the height h and the time t as h' = A^/^/i and t' = \^/^t we can rewrite 
the KPZ equation as 

dh'{-x,t') 



dt' 



= v'y^h'i^, t') + (V/i'(x, t')f + 77(x, t') , 



where i/' = i//(DA^)^/'^ and the noise is still Gaussian with the same strength D, we only consider the 
same two types of perturbations as for the linear systems. 
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We can compare the scaling forms (1531) and (1541) with themselves as well as with 
the scaling forms of the local quantities. We remark that we have: a = —2/3 and 
a£ = —1/3, 6 = — 2 and hi = —2/3, \c = \x ~ ^/^ cind Xc^ = Xxe — 2/3, i.e. the global 
exponents always differ from their corresponding local exponents. 

Finally, let us mention that also for the KPZ case one does not obtain non-trivial 
limits for global fluctuation-dissipation ratios. Still, one can perform a similar analysis 
as done in the previous section for the linear systems, thereby recovering the same three 
regimes, but with a different dependence on t, X ~ t~^/^, in the correlated regime. 

5. Discussion and conclusion 

Our study of global quantities in growth processes has revealed some interesting results, 
especially in the correlated regime, and raises a range of open questions that warrant a 
more in-depth study in the future. 

A prediction of the scaling behaviour of our global quantities is far from obvious, due 
to the complicated nature of the mean square width that involves all Fourier components 
(with the exception of the zero mode) of the Fourier series representation of the surface 
height. We find the available results to be compatible with the scaling form 

C{t,s) = s'^/^+''/^fcit/s), (55) 

with fc{y) ~ y2a/2-i £qj, y ^ 1, as is readily checked by recalling that for the linear 
equations a = ^-^ and z = m, whereas for the one-dimensional KPZ equation a = | 
and z = |. Similarly, for the the global response of the square of the surface width 
to a change in the noise strength (we need to specify this for the linear systems) is 
compatible with the scaling form 

X{t,s) = s'^/^f^{t/s), (56) 

with fxiy) ~ y2a/z-i fQj^ y ':^ I. We are not able to derive these scaling forms at this 
stage. 

In fact, the situation for the response function is rather complex. The reason is of 
course that for the linear systems the global response depends on how the system was 
perturbed, as the different protocols lead to different scaling forms. This is different 
for the non-linear KPZ equation, as here the same scaling forms are found for the 
different perturbations. Obviously, the non-linearity provides an efficient smoothening 
mechanism that allows the system to relax to the state of the unperturbed system in 
a way that is independent on how the system has been perturbed. This mechanism is 
absent in the linear systems, yielding the observed dependence of the response on the 
protocol. A more quantitative study of this effect is left for the future. 

Our study also reveals that in the correlated regime the usual limit value of the 
fluctuation-dissipation ratio involving global quantities only leads to trivial results. 
However, the ratio ( 15T]) at flnite observation times still allows us to distinguish between 
the different growth regimes. 
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In order to measure a global response one needs to perturb in a uniform way the 
system. On the level of the Langevin equations, this type of perturbation can only be 
achieved by varying one of the parameters entering these equations. As we argued in 
Section 2.3, based on the microscopic model, the surface tension u is directly related 
to the quantity conjugated to the square of the surface width. This is different for the 
strength D of the noise correlator. It is therefore interesting to note the similarities in 
the response of the squared width due to a change of either the surface tension or the 
noise strength. In both cases the non-equilibrium exponents that govern the behaviour 
in the ageing regime are not identical to the corresponding exponents describing the 
correlation function of the squared width. This seems therefore to be a general property 
of global two-time quantities. 

Experimentally a global response is much easier to measure than a local one. Local 
perturbations and local measurements are not easily achieved in a growing system. This 
is different for a global response that is readily observed when changing the experimental 
conditions, for example through a change in temperature, as this automatically affects 
the whole system. Another advantage of studying global quantities in the context of 
growth processes is seen at the upper critical dimension. Indeed, global quantities 
continue to display a standard dynamical scaling behaviour at the critical dimension. 
Local quantities, however, generically display a logarithmic dependence on time in that 
case which makes an analysis much more demanding. For these reasons we hope that our 
study will motivate other groups, especially experimental ones, to use global quantities 
for the characterisation of non-equilibrium growth processes. 

Acknowledgments 

We thank Sebastian Bustingorry for bringing Reference [H] to our attention and Uwe 
Tauber for an enlightening discussion on conjugated quantities within a field-theoretical 
setting. This work was supported by the US National Science Foundation through grant 
DMR-0904999. 

6. References 



[1 
[2 
[3 

[4; 

[5 
[6 

[7; 

[8 

[9 

[10 

[11 



Meakin P 1993 Phys. Rep. 235 189 

Halpin-Healy T and Zhang Y-C 1995 Phys. Rep. 254 215 

Barabasi A-L and Stanley H E 1995 Fractal Concepts in Surface Growth (Cambridge: Cambridge 

University Press). 
Krug J 1997 Adv. Phys. 46 139 
Krug J 1995 in Scale Invariance, Interfaces, and Non-Equilibrium Dynamics, edited by A. 

McKane et al. (New York: Plenum, New York) 
Family F and Vicsek T 1985 J. Phys. A 18 L75 
Chou Y-L and Pleimling M 2009 Phys. Rev. E 79 051605 
Kallabis H and Krug J 1999 Europhys. Lett. 45 20 

Rothlein A, Baumann F and Pleimling M 2006 Phys. Rev. E 74 061604 
Rothlein A, Baumann F and Pleimling M 2007 Phys. Rev. E 76 019901(E) 
Bustingorry S, Cugliandolo L F and Iguain J L 2007 J. Stat. Mech. P09008 



Characterisation of non- equilibrium growth through global two-time quantities 24 



[12: 
[13 

[14 

[is: 

[16 

[17: 
[is: 

[19 

[20" 
[21 
[22 
[23 

[24: 
[25 

[26: 

[2?: 
[28: 

[29 
[30 
[31 

[32 
[33 

[34 

[35 

[36: 

[3?: 
[38: 

[39 

[40 
[41 
[42 
[43: 
[44: 
[45 
[46: 
[4?: 
[48: 
[49: 
[50 
[51 
[52 
[53: 
[54 
[55 
[56 

[57: 



Bustingorry S 2007 J. Stat. Mech. P 10002 

Chou Y-L, Pleimling, M, and Zia R K P 2009 Phys. Rev. E 80 061602 

Nguyen T T T, Bonamy D, Phan Van L, Cousty J, and Barbier L 2009 larXiv:0904.2200l 

Henkel M, Noh J D, and Pleimling M 2010 unpublished 

Henkel M and Pleimling M 2010 Non-equilibrium phase transitions Volume 2: Ageing and 

dynamical scaling far from equilibrium (Heidelberg: Springer) 
Mayer P, Berthier L, Garrahan J P, and Sollich P 2003 Phys. Rev. E 68 016116 
Pleimling M and Gambassi A 2005 Phys. Rev. B 71 180401(R) 
Annibale A and Sollich P 2006 J. Phys A: Math. Gen. 39 2853 
Calabrese P, Gambassi A, and Krzakala F 2006 J. Stat. Mech. P06016 
Dutta S B, Henkel M, and Park H 2009 J. Stat. Mech. P03023 
Edwards S F and Wilkinson D R 1982 Proc. R. Soc. London Ser. A 381 17 
Mullins W W 1963 in Metal surfaces: Structure, energetics, and kinetics (Metals Park, Ohio: 

Am. Soc. Metal) 
Kardar M, Parisi G, and Zhang Y-C 1986 Phys. Rev. Lett. 56 889 
Family F 1986 J. Phys. A 19 L441 

das Sarma S, Lanczycki C J, Kotlyar R, and Ghaisas S C 1986 Phys. Rev. E 53 359 
Meng B and Weinberg W H 1996 Surface Science 364 151 
Barteh M C and Evans J W 1999 Surface Science 423 189 
Liu Z-J and Shen Y G 2005 Surface Science 595 20 
Elsholz F, Scholl E and Rosenfeld A 2004 Appl. Phys. Lett. 84 4167 
Ghaisas S V 2001 Phys. Rev. E 63 062601 

Barbier L, Masson L, Cousty J, and Salanon B 1996 Surface Science 345 197 
Dougherty D B, Lyubinetsky I, Williams E D, Constantin M, Dasgupta C, and Das Sarma S 

2002 Phys. Rev. Lett. 89 136102 
Dougherty D B, Tao C, Bondarchuk O, Cullen W G, Wilhams E D, Constantin M, Dasgupta C, 

and Das Sarma S 2005 Phys. Rev. E 71 021602 
Bondarchuk O, Dougherty D B, Degawa M, Williams E D, Constantin M, Dasgupta C, and Das 

Sarma S 2005 Phys. Rev. B 71 045426 
Horowitz C M, Monetti R A and Albano E V Phys. Rev. E 63 066132 
Chame A and AarSo Reis FDA 2002 Phys. Rev. E 66 051104 
Horowitz C M and Albano E V 2003 Eur. Phys. J. B 31 563 
Kolakowska A, Novotny M A and Verma P S 2004 Phys. Rev. E 70 051602 
Muraca D, Braunstein L A and Buceta B C 2004 Phys. Rev. E 69 065103(R) 
Braunstein L A and Lam C-H 2005 Phys. Rev. E 72 026128 
Irurzun I, Horowitz C M and Albano E V 2005 Phys. Rev. E 72, 036116 
Horowitz C M and Albano E V 2006 Phys. Rev. E 73 031111 
Aarao Reis FDA 2006 Phys. Rev. E 73 021605 

Oliveira T J, Dechoum K, Redinz J A and Aarao Reis FDA 2006 Phys. Rev. E 74 011604 
Kolakowska A, Novotny M A and Verma P S 2006 Phys. Rev. E 73, 011603 
Majaniemi S, Ala-Nissila T, and Krug J 1996 Phys. Rev. B 53 8071 
Giesen M 2001 Prog. Surf Sci. 68 1 

Dougherty D B, Lyubinetsky I, Einstein T L and WiUiams E D 2004 Phys. Rev. B 70 235422 
Wolf D E and Villain J 1990 Europhys. Lett. 13 389 
Golubovic L and Bruinsma R 1991 Phys. Rev. Lett. 66 321 
Das Sarma S and Tamborenea P 1991 Phys. Rev. Lett. 66 325 

Iguain J L, Bustingorry S, Kolton A B and Cughandolo L F 2009 Phys. Rev. B 80 094201 
Crisanti A and Ritort F 2003 J. Phys A: Math. Gen. 36 R181 
Cughandolo L F, Kurchan J, and Peliti L 1997 Phys. Rev. E 55 3898 
Fielding S and Sollich P 2002 Phys. Rev. Lett. 88 050603 
Calabrese P and Gambassi A 2004 J. Stat. Mech. P07013 



Characterisation of non- equilibrium growth through global two-time quantities 25 



[58 
[59 
[60 
[61 
[62 
[63 
[64 
[65 



Mayer P, Berthicr L, Garrahan J P and SoUich P 2006 Phys. Rev. Lett. 96 030602 

Garriga A, Pagonabarraga I, and Ritort F 2009 Phys. Rev. E 79 041122 

Newman T J and Bray A J 1996 J. Phys. A: Math. Gen 29 7917 

Lam C-H and Shin F G 1998 Phys. Rev. E 58 5592 

Buceta R C 2005 Phys. Rev. i; 72 017701 

Miranda V G and Aaro Reis FDA 2008 Phys. Rev. E 77 031134 

Wio H S, Revelli J A, Deza R R, Escudero C and de la Lama M S 2010 EPL 89 40008 

Newman T J and Swift M R 1997 Phys. Rev. Lett. 79 2261 



